Quantum-classical correspondence of the Dirac matrices: 
The Dirac Lagrangian as a Total Derivative 
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The Dirac equation provides a description of spin 1/2 particles, consistent with both the principles 
of quantum mechanics and of special relativity. Often its presentation to students is based on 
mathematical propositions that may hide the physical meaning of its contents. Here we show that 
Dirac spinors provide the quantum description of two unit classical vectors: one whose components 
are the speed of an elementary particle and the rate of change of its proper time and a second vector 
which fixes the velocity direction. In this context both the spin degree of freedom and antiparticles 
can be understood from the rotation symmetry of these unit vectors. Within this approach the 
Dirac Lagrangian acquires a direct physical meaning as the quantum operator describing the total 
t ime-derivat ive . 
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I. INTRODUCTION 



Historically, Paul Dirac found the Klein-Gordon equation physically unsatisfactory for the appearence of nega- 
tive probabilities arising from the second-order time derivative [ij. Thus he sought for a relativistically invariant 
Schrodinger-like wave equation of first order in time derivative of the form (h = c = 1): 

idt^p^H^. (1) 

In order to have a more symmetric relativistic wave equation in the 4-momentum components, he also looked for an 
equation linear in space derivatives i.e. momentum P = — «V, so that H takes the form 

H = a-P + Pm, (2) 

with a and /? being independent on space, time and 4-momentum. The condition that Eq. ^ provides the correct 
relativistic relationship involving energy, rest mass and momentum, 

= p^ + , (3) 

requires that a and /? = 0:4 obey the anticommutation rules {ai,aj} — 2Sij (« = 1,4) 0]. Dirac found that a set of 
4x4 matrices satisfying this relation provides the lowest order representation of the four a^. They can be expressed 
as a tensor product of (2 x 2) Pauli matrices pi and ai belonging to two different Hilbert spaces: ai = pidi (i = 1,3) 
and (3 = 03 l|- As a consequence tp in Eq. ([1]) is a 4-component wave function. It turns out that a is the velocity 
operator [ij, [3 1 . So that tp''aip is the current density, determining the coupling with the electromagnetic field [1, 2, ^ 

This derivation, although straightforward, does not help for a clear physical understanding of the four component 
wave function and of Dirac matrices. Indeed Richard Feynman in his Nobel lecture pointed out that Dirac obtained 
his equation for the description of the electron by an almost purely mathematical proposition. A simple physical view 
by which all the contents of this equation can be seen is still lacking. 

The fact that relativistic Dirac theory automatically includes spin leads to the conclusion that spin is a purely 
quantum relativistic effect originating from the finite dimensional representations of the Lorentz group Never- 
theless this interpretation is not generally accepted, e.g. following Weinberg argument (Ref. [3| Chapter 1) . . .it is 
difficult to agree that there is anything fundamentally wrong with the relativistic equation for zero spin that forced 
the development of the Dirac equation - the problem simply is that the electron happens to have spin h/2, not zero. 
Technically speaking, the homogeneous Lorentz group (in contrast e.g. to the group of rotations) is not a compact 
group, hence the implementation of this symmetry in quantum mechanics does not need the use of finite dimensional 
representations. On the contrary, being non compact, it has no faithful finite dimensional representation that is 
unitary [J, Thus the homogeneous Lorentz group is the only group of relativistic quantum field theory acting 
on multiple-components quantum fields non-unitarily. This rather surprising fact conflicts with a theorem proved 
by Wigner in 1931 (see Ref. [2j Chapter 2) which states that any symmetry operation on quantum states must be 
induced by a imitary (or anti-unitary) transformation. 
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The conflict is usually overcome, either by regarding the field not as a multicomponent quantum wavefunction but 
as a classical field (JL or by pointing out that the fundamental group is not the (homogeneous) Lorentz group but 
the Poincare group Independently on the point of view, one consequence is that the Hermitean conjugate -0^ of 
the (four-component) spinor field ip does not have the inverse transformation property of as requested by quantum 
mechanics. The rather ad hoc, though generally accepted, solution is to define — -0^7° called the Dirac conjugate 
of 0, being 7° = /? the time Dirac matrix 0, 0, H 01 • 

In this paper we address some naturally arising questions: is there a fundamental compact symmetry group that 
requires the occurrence of spin? is there a reason why the velocity operator a of elementary matter particles is a 
vector made of 4 x 4 matrices instead of continuous variable operators? is there a simple interpretation of Dirac 
four-component spinors? Here we discuss how the spin of elementary particles can be understood as a consequence of 
a rotation symmetry displayed by the kinematic variables describing their space-time motion. Furthermore we show 
that such rotation symmetry implies the existence of antiparticles. 

Before starting our analysis, it is worth pointing out that here we follows the Dirac old fashioned point of view, 
regarding the Dirac equation as a Schroedinger-like quantum mechanical equation providing the quantum description 
of a relativistic point-like particle. In contrast the modern view regards the Dirac equation as a classical wave equation 
describing the 1/2 spin field. The field is then quantized by means of canonical quantization for relativistic fields 



II. ROTATIONS IN QUANTUM MECHANICS 



A rotation in classical physics is implemented by a 3 x 3 orthogonal matrix R which, acting on a given vector v, 
gives the rotated vector v' = i?v. Rotations are a symmetry group whose generator is the angular momentum [§]. In 
quantum mechanics a symmetry transformation is a transformation of the state kets describing the physical system 
and of the operators O in the ket space. The transformation operators D{R) are unitary and have the same group 
properties as R (see also supplementary information). The expectation values of the angular momentum operators 
transform under rotation as classical vectors: 

(Jfe) ^ Rki (Ji) (4) 

where repeated indices are implicitly summed over . The lowest number iV = 2j -I- 1 of dimensions in which the 
angular momentum commutation relations are realized is = 2 (j = 1/2). In this case the angular momentum 
operators can be represented in terms of the Pauli matrices: Jk = crk/2. Independently on the physical state in the 
2D Hilbert space, they obey the following relationship: 

E(^')' = 1 (5) 

i 

which is not satisfied by higher order angular momentum operators. Hence Pauli matrices are the best quantum 
correspondents of classical unit vectors u (with — 1) [10]- there is any classical physical variable which 

is described by a unit vector, Pauli matrices provide its natural quantization: Uj — > cTj, which preserves rotation 
symmetry and ensures expectation values which maps on the classical values. 



III. FROM RELATIVISTIC KINEMATICS TO UNIT VECTORS 



Let us consider a relativistic point-like particle moving with velocity v{t) with respect to an observer's reference 
frame K 11]. A possible spacetime worldline trajectory is shown in Fig. la. We can consider the Lorentz- invariant 
quantity 

dr(t)2 = fit2 _ ^x(i)2 . (6) 

In the instantaneous rest frame K' of the particle dx = and hence dT{t)^ = dt^ which admit solutions dr = ±dt. 
The solution dr = dt is the one usually presented in the introductory texts on special relativity (see e.g. [HI). More 
recently Costella et al. [l^ pointed out that the negative solution provide a classical description of antiparticles in 
agreement with the Feynman's interpretation of antiparticles as particles going hack in time [7[. 
Thus from Eq. ^ it follows that: 



(7) 
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where 7 = (1 — w^)"^/^ is the usual Lorentz boost parameter yielding time dilation. 

The time r is the proper time of the particle. Thus the quantity dr/dt = t = zLl/j(t) expresses the rate of change 
of the proper time with respect to the time t of the reference frame i.e. proper time velocity. 

The usual way to describe spacetime trajectories (worldlines) relies on i — x Minkowski diagrams, as shown in fig. 
Illlb ,. where the light cones separates causal events from space-like events. Additional information can be inferred by 
drawing proper time-time trajectories (timelines) as shown in fig. Illlb . The solid line describes the trajectory of a 
particle with f > 0. Trajectories with r > can be identified by classical particle motion states, while solutions with 
f < correspond to classical anti-particle states . 




FIG. 1: Worldlines of relativistic point-like particle and antiparticles. a) In the usual Minkowski space-time diagram particles 
and antiparticles worldlines are not distinguished, b) A proper time-time diagram manifests the different kinematics of particles 
(r > 0) and antiparticles (f < 0) with respect to an observer's reference frame. 



When attempting to understand the spin and Dirac spinors on a physical ground, the following question arises: are 
there kinematic variables which can be described as components of unit vectors ? 

The speed of a particle v{t) — |x(f)| is limited (by the speed of light c — 1) like the component of a unit vector. 
Additionally the quantity f can be viewed as the complementary component of this unit vector. Indeed by definition 
of 7, it follows that: 

Ixl' + r^-l. (8) 

Thus we can introduce a convenient miit vector f = (|x| , 0, t) lying on the rir^ plane (see Fig. Illlb ,) and an angle (f) 
so that: 

|x|=sin(/), T — cos(j). (9) 

A particle changing its kinematic state will move on the unitary circle defined by the angle 0. Figure Ullk displays 
one such kinematic vector with positive components. If the particle is at rest with respect to the reference frame, 
the vector lies on the r^ axis (f = 1). Moreover f decreases when the particle-speed increases as predicted by special 
relativity (time dilation). 

The unit vector f, besides f, is able only to describe the modulus of the particle velocity v — |x|. The particle 
velocity is indeed a 3D vector and the direction of x can be accounted for by one additional 3D unit vector s providing 
just the direction of x. Hence the motion state of a particle can be described by a specific couple of unit vectors f 
and s: 

T==r3, x = ris. (10) 

Within this approach changes of the particle velocity respect to an inertial frame can be accounted for by rotations 
of f in the kinematic plane (changes of the modulus) and rotations of s (changes of the direction). The unit vector s 
can be transformed according to arbitrary 3D rotations around an arbitrary 3D unit vector n: 



Si s'i = [RniO)]ij Sj, 



(11) 
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FIG. 2: Geometric representations of unit vectors defining the kinematic state of a relativistic point-like particle, (a) The unit 
vector f is defined by the modulus of space and proper time velocities. A particle changing its kinematic state will move on 
the unitary circle defined by the angle (j). (b) The unit vector s is defined by the spatial direction of velocity and can rotate on 
the 3D unitary Poincare sphere. 



where 6 labels the angle of rotation about n. On the other hand, physical kinematic states f admit only rotations 
about the r2-axis: 

n^r',^[B!^m,,r,. (12) 

We observe that this rotation symmetry suggests that states obtained rotating f , and § should be considered 
as accessible states. In particular f also describes kinematic states in the second and third quadrant with f < 0, 
corresponding to classical antiparticles [l3| 

From the point of view of classical special relativity, the description in terms of s and f appears to be redundant: a 
given velocity is described by two different states. For example a particle with given velocity along direction d = x/ |x| 
can be described by the unit vectors S| = d and f| = (sinc^, O,cos(^) with </> = arcsin|x|, or equivalently by the unit 
vectors S| = — d and f ^ — (sin </)', 0, cos (j>') with 0' = —(f). This two-valuedness can be described in terms of the helicity 
variable h = (s-p)/p. Below we show that this classical twofold degeneracy is indeed the classical correspondent of the 
helicity states determined by quantum spin. It is worth pointing out that, although the present approach describes a 
spin-like degree of freedom yet at a classical level, the interaction of a classical particle with the electromagnetic field 
is not affected by this additional degree of freedom, in contrast to what happens after quantization. 

Figure 3 provides a clear geometric interpretation of the different kind of kinematic states: the first quadrant 
contains spin up particles, the second one spin up antiparticles, the third spin down antiparticles and the fourth 
quadrant spin down particles (see Supplementary Information). 



IV. FROM UNIT VECTORS TO DIRAC MATRICES 



As discussed earlier, a quantum mechanical description of unit vectors is obtained replacing the classical vector 
components with Pauli matrices: — > pi and Si at, where pi and at are now Pauli matrices acting on two different 
Hilbert spaces TZ and S. Hence: 

T = r-i ^ p3 ^ I = (3 

X = riS^pi(g)cr = Q;, (13) 

being / the identity operator in the s space describing the direction of velocity. The a and /? are indeed the well-known 
Dirac matrices in the standard representation Thus we find a direct symmetry-driven correspondence showing 
that a is the velocity vector operator, and P the proper-time velocity operator. In this quantum picture the helicity 
operator for a particle of given momentum becomes /i = (c • p)/p in agreement with the Dirac theory [J and with our 
classical derivation of helicity states. In summary, at a classical level, the rotation symmetries wc introduced describe 
the correct relativistic kinematics. Furthermore they lead to two possible helicities for a given motion state and the 
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FIG. 3: Sections of the velocity space describing the different solutions. Unit vectors in the first and fourth quadrants describe 
spin up and spin down particles respectively. Unit vectors in the second and third quadrant describe spin up and spin down 
antiparticles. 



existence of particles with f < 0. After quantization, without any other additional assumption, they give rise to the 
4-component Dirac spinors fully describing quantum spin and antiparticles with the corresponding expectation value 
{(3) < 0. We point out that in this framework the Dirac matrices are derived only on the basis of the fundamental 
rotation symmetry without any reference to the Dirac equation or to the energy of a relativistic free particle. 

At this point it is interesting to show that the energy of a relativistic free particle can be written as observed By 
Breit as -E = v ■ p + rrn. This equation has the same structure of the Dirac Hamiltonian. 



V. THE DIRAC LAGRANGIAN AS A TOTAL DERIVATIVE 



A crucial feature of this symmetyry based derivation of the Dirac matrices is that it has been obtained by regarding 
the position-coordinate and the proper time of the particle on the same footing as functions of the time-coordinate 
of the reference frame, i.e. we use t as the relevant meter to which compare the dynamical evolution of all other 
observables (as happens in non-relativistic mechanics). Accordingly if a quantum wavefunction depends on position 
X, it should also depend on r i.e. 'ip{t,x.,T). We now consider a particle with such wavefunction and require the 
conservation of its probability density P(i, x, r) = tp*(tj ''")V'(ij ''")■ 

Stationarity of ip (and ip*) ensures that a continuity equation holds. In turn if, as in Lagrangian formulations with 
complex fields, ip and ip* are regarded as independent, Eq. p4p implies stationarity of ip (and of ip*) i.e. dip / dt = 0. 
Because the time dependence of the wavefunction is both explicit and implicit (through the time dependence of x and 
r), we have: 

^^ = dtTp + ±-d^lp + Tdrlp =0. (15) 

dt 

Thus the request of stationarity and the crucial symmetry-based quantization Eq. ([TH)) . provides a generalized Dirac- 
like equation: 

idtipit, x, t) = Hg%l){t, X, t) , (16) 

with 

Hg = a- (-z^x) + P{-idr) . (17) 
The corresponding Lagrangian density can be written as 

ip^Cip = ijh{dt + a-d^ + (3dr)tlj . (18) 
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This Lagrangian operator C has a precise physical meaning, being the Hermitean quantum operator describing the 
total time derivative: id/dt C after symmetry-based quantization Eq. p^ . Equation (|17p is more general than 
the Dirac equation since the mass parameter is now replaced by the internal time momentum operator. Elementary 
matter particles of given mass can be viewed as eigenstates of this operator. In this case Eq. pT|) reduces to the 
standard Dirac equation. On the other hand, this new degree of freedom offers a chance to shed new light on some 
fundamental aspects and concepts. 

As an example this new degree of freedom could be exploited to discuss charge conservation. Charge conservation 
is generally accounted by invariance under phase change of a field, Eq. (|17p allows the interpretation of this phase 
change as the consequence of proper time translation applied to an eigenstate of the proper time operator. Charge 
conservation could thus be viewed as arising from invariance of the Lagrangian under proper time translations. 

It is also interesting to address the proper time reversal symmetry r — r. Applying this symmetry to a positive 
energy solution of Eq. [T7| which is also an eigenstate of the proper time momentum (i.e. a particle with fixed mass 
and charge), we obtain its antiparticle with positive energy and opposite charge, so antiparticle solutions with positive 
energy emerge without the need for second quantization adjustments. 

VI. CONCLUSIONS 

Within the approach here presented, changes of the velocity modulus and direction of a relativistic particle can 
simply be accounted by rotations of two independent unit vectors. Dirac spinors just provide the quantum description 
of these rotations. These transformations are able to describe the helicity of a pointlike particle yet at a classical 
relativistic level, rendering spin a less mysterious degree of freedom. This analysis sacrifies explicit covariance by 
making explicit rotation symmetry which nevertheless is the fundamental symmetry on which the algebra of the 
Lorentz group is based. We observe that the addressed rotation symmetries form a compact group hence with unitary 
finite-dimensional representations as all other symmetry groups in quantum field theory. A feature of this analysis is 
that it has been obtained by regarding the position-coordinate and the proper time of the particle on the same footing 
as functions of the time-coordinate of the reference frame, i.e. t is used as the relevant meter to which compare the 
dynamical evolution of all other observables (as happens in non- relativistic quantum mechanics). 

Within this approach we derived the Dirac equation just invoking total stationarity of the wavefunction with 
respect to the reference time. No assumptions about the classical relativistic energy of a particle or about the 
quantum operator replacements {E —> idt and p — j(?x) have been performed. Finally we observe that the quantum 
replacement m — > —idr implies an internal time-energy uncertainty principle At Am, e.g. as required by a gedanken 
experiment recently proposed by Aharonov and Rezni p^ . 
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